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ABSTRACT 

Pointwise  bounds  of  a bivariational  nature  are  derived  on  the  solution  of  a standard 
Fredholm  integral  equation  of  the  second  kind,  with  a symmetric  kernel.  The  bounding 
functionals  involve  two  trial  vectors,  one  approximating  the  solution  and  the  other  ap- 
proximating the  reciprocal  kernel.  Even  with  the  latter  taken  as  the  null  vector,  there 
is  significant  improvement  over  a previous  approach.  It  is  shown  how  suitable  choices  of 
trial  vector  lead  to  expressions  for  bounds  on  Neumann,  Pade  and  Fredholm  approximations. 
AMS(MOS)  Subject  Classification  - 45L05,  45G20 

Key  Words  - Integral  equation.  Symmetric  kernel,  Pointwise  bivariational  bounds, 

Hilbert  space,  Neumann,  Pade,  Fredholm. 

Work  Unit  Number  1 - Applied  Analysis 

EXPLANATION 

It  is  not  often  possible  to  solve  exactly  an  integral  equation  like 

b 

<|>  (x)  + A/  k (x,y)  (y)dy  = f (x)  , a <_  x £ b , 

for  the  unknown  function  <t>(x)  in  terms  of  the  known  functions  f(x)  and  k(x,y)  (the 
'kernel'  function,  symmetric  in  x and  y) . So  we  derive  analytical  expressions  for  two 
other  functions  4>_  ( x)  and  <J>+  ( x)  with  the  property 

4>_(x')  £4>(x')  £<f+(x'),  a <_  x'  £b  , 

providing  upper  and  lower  bounds  on  ( x ' ) at  all  points  of  the  interval.  The  differences 

<}r(x')  - <}>  ( x ' ) and  <f>(x')  - 4>_  ( x ' ) have  the  same  kind  of  size  as  A 6(x'),  where 
t On  leave  from  Bradford  University,  England. 
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POINTWISE  B I VARIATIONAL  BOUNDS  ON  SOLUTIONS  OF  FREDHOLM  INTECRAL  EQUATIONS 

Peter  D.  Robinson^ 


1.  Introduction. 

In  this  paper,  pointwise  bounds  of  a bivariational  nature  are  derived  on  the  solutions 
of  a standard  Fredholm  integral  equation  of  the  second  kind,  with  a symmetric  kernel.  The 
method  depends  on  improved  bivariational  bounds  associated  with  linear  equations  in  a Hilbert 
space,  and  exploits  the  similarity  between  the  integral  equation  itself  and  the  equation 
specifying  the  reciprocal  kernel.  The  bounding  functionals  involve  two  trial  vectors,  one- 
approximating  the  solution  of  the  integral,  equation  and  the  other  approximating  the  reciprocal 
kernel.  Even  when  the  latter  is  taken  to  be  the  null  vector,  the  bounds  refine  those  which 
can  be  obtained  from  an  approach  of  Lonseth  (1);  a simple  numerical  illustration  is  given.  It 
is  also  shown  how  suitable  choices  of  trial  vector  lead  to  expressions  for  bounds  on  approxi- 
mate solutions  of  Neumann,  Pad£  and  Fredholm  type. 

2 . Hilbert  Space  Formulation. 

The  pointwise  bounds  are  obtained  as  instances  of  complementary  (upper  and  lover)  bi- 
varational  bounds  on  the  linear  product  <g,$>  associated  with  a pair  of  equations 


H - f , 

Alp  - g , 


ip , f e 3f 
ip,g  e 3( 


in  a Hilbert  space  7(  , A being  a self-adjoint  operator  on  3C.  We  assume  that  A is  strict)'- 
positive  and  bounded  below  away  from  zero,  so  that  for  some  positive  number  8 

«>,A$>  >.  8 <4’,4,>  for  all  4>  e JC  . (2.3) 

Then  the  inverse  operator  A ^ exists  with  domain  the  whole  of  Jf.  Further,  we  assume-  that 
A is  bounded  above,  so  that  for  some  positive  number  a 

>_  <■  t^A^^for  all  t c JC , a > 8 > 0 • (2.4) 

(The  case  a -*■  m is  admissible,  but  means  that  A is  not  bounded  above;  then  A is  merely 
defined  in  JC  and  wherever  necessary  vectors  must  be  assumed  to  belong  to  the  domain  of  A)  . 

The  Fredholm  integral  equation  with  real  symmetric  kernel 

♦ (x)  + \ /^(x.yHfyJdy  = f(x),  a < * < b , (2.5) 
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is  identified  with  (2.1).  If  i|i(x,y)  is  the  reciprocal  (or  resolvent)  kernel,  so  that  the 
relation 

b 

♦ (x)  + A / iMx.y)  f(y)dy  = f(x) , a <_  x < b , (2.6) 

specifies  the  solution  $(x)  of  (2.5),  then  it  is  a standard  result  that 
b 

<’(x,x,)+A/ak(x,y)iHy,x,)dy  = k(x,x'),  a£x£b,  a<x'<b.  (2.7) 

This  equation  (2.7)  is  identified  with  (2.2).  We  take  for  Jf  the  real  space  of  square- 
inteqrable  functions  with  inner  product 

b 

<hl,h2>  = /a  h1<x^  h2(x)dx,  for  all  h ,h2  e ^ ' (2.8) 

although  it  is  possible  to  adapt  the  theory  for  complex  spaces. 

The  symmetric  kernel  k(x,y)  is  assumed  to  be  of  Hi lbert-Schmidt  type,  and  describes 
a compact,  self-adjoint  operator  K in  3C,  so  that 

A = I + AK  . (2.9) 

The  reciprocal  kernel  iMx,y)  is  also  symmetric,  and  describes  the  reciprocal  operator 
T given  by 

r = K(I  + AK)”1  = KA_1  . (2.10) 


Accordingly  complementary  bounds  on  < g , <J>  > lead  to  complementary  bounds  on  <Mx'),  for 
any  x*  in  a ± x*  <_  b . 

The  strict  positivity  requirement  (2.3)  places  a restriction  on  the  values  of  the 
parameter  A for  which  the  results  will  be  valid.  Since  K is  compact  and  self-adjoint, 
it  follows  that  for  some  non-negative  numbers  M and  L, 

M < *,♦  > ^<$,K*>  >_  -!.<♦,*>  for  all  * c Jf  . (2.13) 


-2- 


In  theory,  the  numbers  +M  and  -L  are  respectively  the  positive  and  negative  eigen- 
values of  K of  greatest  magnitude,  but  in  practice  they  can  be  taken  as  the  best  avail- 
able upper  and  lower  bounds  to  these  eigenvalues.  Then  we  assume  that 


L 1 > X > -M_1 


(2. 1C) 


and  take 


a = 1 + XM,  6 = 1 - XL  if  X > 0 , 
a = 1 + XL,  (3  = 1 + AM  if  A<0. 


(2.15) 


The  bi variational  bounds  derived  in  §3  below  involve  the  positive  constants  Q 


and  n defined  by 


K = j(B  1tu  1),  n = j(6  X-a  X)  . 


(2. 1C) 


For  either  positive  or  negative  X,  it  follows  from  (2.15)  that 


1 + |X(M-L) 

f'  ~ "(i+XMr (1-xl)  ' 


j|x| (M+L) 

(I+xmhT^xl)  ‘ 


(2.17) 


3 . Comp 1 erne nl . ary  Biva riation.il  Bound s 
Associated  with  an  equation 


AO  = h,  0,h  c If  , 

who>;e  solution  0 is  unknown,  are  the  complementary  variational  bounds  on  <h,fi  > : 

GQ((-);h)  < <h,0>  £ Gg((-',h),  for  all  O £ If  , I 


G (C-);  h)  = -<  0,AfJ>  + 2<  h,0>  + a_1  ||aO  -h  | 


G(,(0;h)  ■-=  -<0,AC-)>  +2<h,0>  -»  B-1  |(aC-)-1, |(  2 . (3-',) 

The  norm  here  is  the  usual  Hilbert  space  norm,  ||  <>||2  = < . t)io  variational  boun' 

ing  properties  of  g and  G follow  at  once  from  the  identities 
a p 


G (W;h)  — < h , 0 > - a_1<A1/?60,  (a-A)  A1/2fi0>  , 


(3.5) 


and 


G (U;h)  = < h , 0 > 4 p”1  <A)/2  60,  (A-p)A1/260>  , 

in  terms  of  the  difference-vector  60  = 0 - 0,  together  with  (2.3)  and  (2.4). 
Consider  now  the  pair  of  equations 

A(s4>  + ti(i)  = sf  4 tg  , 


A(s 6 - tiJO  = sf  - tg  , 


obtained  from  (2.1)  and  (2.2),  s and  t being  scalar  multipliers.  If  we  subtract  the 
two  inner  products  like  <h,  0>  which  can  be  bounded  as  in  (3.2)  we  obtain: 

<(sf4tg)  , (s<j>  4 tip)  > - < (sf-tg  ) , (s<f>-t^)  > = 2st<g,  <j>>  4 2st<f,iJ> 

= 4st<g,<J>>  • (3.9) 

(The  last  step  follows  since  <f,t|‘>  = < A <> , ti-'  > = = <$,g>  - <q,ij.>)  . Thus, 

from  (3.2)  and  (3.9),  we  have  for  all  <!> , H*  c If  : 

(sf  4 tg;  s4>  4 tV)  - G^,  (sf  - tg;  s4>  - tV)  £ 


< 4st<q,  < G„  (sf4tg ; s4  4t'i')  - G (sf-tg;  si'-tV) 

— — b a 


(3.10) 


The  lert-hand  and  right-hand  members  of  (3.10)  reduce  to 


4st[  <T,A<f'>  4 < 4'  ,f>  4 <g,4>}  4 


4£,st<A4-f  ,A¥- g > + 2n{s2  | j A4>—  f ||  4 t.2  ||  AV-g  j|2} 


(3.11) 


Dividing  (3.10)  through  by  4 | st | , and  choosing  the  optimal  ratio 


s | A4-f  ||  = 1 1 1 II  AV-g  ||  , 


(3.12) 


we  obtain  the  fin  d re.  It 


J ( V , <! ) 4 CS(i',<f)  - nC(¥,*)  < <g,4>  < J (¥,«>)  4 (S(¥,4>)  4 l!C(i  ,4') 


(3.13) 


j(¥,4)  =■  -<¥,  Af  > 4 < V , f > 4 < g , 4 > - <g,<(>  - , 


(3.14) 


S (V , t ) -•  <A<-f,  AV-g>  w <A6),A64>  , 


(3.1M 
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C(S',*>  = ||AV-<l||  ||  AV-t  ||  = ||  ASlI'll  ||  A6s'.||  (3. 10) 

and 

■?$  = «■  - <*■  e Jf,  c 3C  . (3.17) 

The  expression  of  the  functionals  J,  S and  C in  terms  of  the  difference-vectors  (3.17) 
makes  clear  the  bi variational  nature  of  the  bounds  in  (3.13).  They  are  tighter  than  others 
previously  obtained  [2J  by  exploiting  the  identity 

< (sf  ± s ^g)  , (s<Ji  ± s *|||)  > - s2<f,i}>>-  s 2<g,iJj>  = ± 2<g,4>>  (3.10) 

rather  tlian  (3.9);  see  also  [3J  ■ 

4 . Application  to  the  Integral  Equations 

When  applied  to  the  pair  of  integral  equations  (2.5)  and  (2.7),  the  functionals  in 
(3. 13) -( 3. 16)  take  the  forms 

J(x’)  = X_1{f(x')  - *(>;'))  - <6if',(  1 3 XK)f<j> 

b b 

= - / M1  ( x , x ' ) { ( 1+XK)  4 (x)  - f (x)  )dx  •+  / k (x , x 1 ) 4(x)dx  , (4.)) 

1 a ‘ a 

S(x')  = < (1+XK)(S3V,  (1+Xk)6^> 

b 

= / { (1+XK)  4(x)  - f (x)  }{  (HXK)  V(x.x')  -U(x,x')}dx  , (4.2) 

1 a 

and 

C(x’)  = II  (1+XK) 6$ ||  II  (1+Xk) 6i/i || x 

b « 1 ^ p W") 

= [/  { ( l+X  K)  4 (x)  - f(x)}  dx)  ' \J  { (1+XK)  (x,x* ) - k(x,x')  } dx]  ' . 

1 a 'a 

(4.3) 

Whenever  there  is  ambiguity,  the  operator  K is  understood  to  act  on  the  first  argument 
of  the  subsequent  function,  so  that 

b 

CT(x,x')  - / k(x,y)  V ( y , x ' ) dy  . (4.4) 

a 

Substituting  from  (2.12)  into  (3.13),  v;e  obtain  the  puintwico  hounds 

(.f(x’>  > (x-)  > «._<x'),  a < x'  <_  b , (4.9) 

where 


« 


♦+(x')  = f(x’)  - X{j(x')  + c S(X*))  + I x |n  C(x')  (4.6) 

and 

»_(x*)  = f(x')  - Xlj(x')  + SS(x')}  - I X I n C (x ' ) . (4.7) 

The  functions  4(x)  and  V(x,x'),  which  determine  J(x'),  S(x')  and  C(x’),  are 
intended  to  be  approximations  to  the  exact  solutions  4>(x)  and  i Hx,x‘)  of  equations  (2.5) 
and  (2.7).  The  upper  bound  <f>+  ( x * ) can  be  minimized  with  respect  to  any  disposable  param- 
eters in  4>  and  T (for  example  multipliers  of  some  set  of  basis  vectors),  and  likewise 

♦ (x‘)  can  be  maximized.  Thus  we  have  a method  for  determining  accurate  bounds  on  tf(x'); 
an  example  is  presented  in  §5. 

A simpler  alternative  procedure  is  to  forego  individual  optimization  of  <}>+  and 
and  write  (4.5)  in  the  form 

0 < |«Mx*)  - tf(x')  - X{j(x')  + £S(x’)}l  < | X |tj  C(x’)  . (4.8) 

This  gives  a pointwise  measure  of  the  accuracy  of  the  function 

<Mx’)  = f(x')  - X{  J (x‘ ) + CS(x')  } (4.9) 

as  an  approximation  to  ( x ’ ) . Disposable  parameters  might  then  be  chosen  to  minimize  the 
st-squarcs"  function  C(x')  . Bounds  on  various  analytical  approximations  to  i,  such 
Neumann,  Pade  and  Fredholm  approximations,  can  also  be  obtained  from  (4.0).  By 
litable  choice  of  the  trial  vectors  t and  V it  is  possible  to  coax  f(x')  -XJ(x')  into 
the  form  of  one  or  other  of  these  approximations.  (See  §6  below) . 

5.  Simple  Trial  Functions 

Sensible  results  are  obtained  from  the  foregoing  theory  even  with  the  uninformed 
choice  of  trial  functions  which  are  identically  zero.  When  both  tlx)  = 0 and  Y(x,x')  «=  0, 
(4.8)  gives  rise  to  the  pointwise  bounds 

0 < |<Mx' ) - [f(x')  - XCKf(x’)]  | £ | A |n  ||  f ||  ||k  (x , x ' )|1  . (5.1) 

With  zero  f(x,x'),  but  arbitrary  <t>(x) , (4.6)  and  (4.7)  become 

♦ +(x')  = f(x')  - XK$  (x  ' ) + XCK{t(x')  + XK*(x')  - f (x') } ± I X In  )|(1+XK)4  - f|| 

||  k(x,x')||x  . (5.2) 
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Louse th  U)  has  in  effect  pointed  out  that  in  any  case  *(x')  lies  between 

f(x')  - XK*(x')  ± | A | ||  6*|!  |i  h(x,x')|lx  ; 

bounds  on  the  Hilbert  space  norm  ||  6<f||  then  lead  to  pointwise  bounds  on  , (x  ) • 

(2.3),  these  would  be 

f (x' ) - XK*  (x ' ) ± | X | 6 1 ||  (1  + XK)  * - f |!  ||  k(x,x'  )]lx 
with  B given  by  (2.15),  which  is  evidently  a coarser  result  than  (5.2)  since  f 
from  (2.16). 

lf  <Hx' ) is  taken  as  cf(x')  in  (5.2),  and  the  factor  c is  optimized  to 
the  best  upper  and  lov.’er  bounds,  those  turn  out  to  bo 

2 2 1/2 

*H(x')  = 4>0  (x ' ) + UjPtx')  i Uj{r  (x* ) - p(x')  } , 


(5.3) 
Using 

(5.4) 

> 2r, 


(5.5) 


whore 


*0 


(x’l  ==  f(x')  - X.^Kf(x')  , 


p(x’)  = C(XK  + XV>  f(x’)  - XUf(x') 

r(x')  = n||  (HXK)f||  ||  Xk||x  , 

_ <f,  (1-vXK)  f> 

1 ||  (1+xk)  r|/ 

2 _ J]fj|f <f,  (HXK)f>2 

||a+),K)f||2  II  (1+XK)  f II 4 


(5.0 


But  in  order  to  export  fully  the  bivar rational  nature  of  the  pointwise  bounds  in  (I.O- 
(4.7),  the  crude  choice  V(x,x‘)  » 0 must  be  relaxed.  If  we  choose 


*(x')  cf(x’),  S'fx.x')  = d): ( x , x 1 ) 


(5.7) 


in  (4.1)- (4.  3),  and  simultaneously  determine  c and  d in  each  case  to  optimize  f.o  he.!.. 
(4.0)  and  (*.7),  those  become 


rfx  =4*  + 

vi  s 0 


2 (p-qn.) (p-qv  )1 

p i i 


2 


, 2 2, 

q(t  -q  ) 


(5.8) 


Ut2-q2)u2  - (p  q>.1)2)1/2  ((t2-q‘!)v2  - (p-qVl)2)1A 

t -q  ' 
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Here,  in  addition  to  (5.6),  we  have  used  the  notation 

q(x*)  = + 2A2K2  + A3K3)f  - (Ak  + a\2>  f(x')  . 


t(x')  = n||(l+XK)  f ||  ||  Ak  +A2k 


2 x ' 


< Ak  , Ak  +A  k > 

v (x')=  5 ~L 

II  *k  tA2k2||  2 


2(  i%  II  *k  ||2  < Ak  , Ak  + A2k2>2 

^ if*  +A2k2||2  " II  Ak  +A2k2||J 


wherein  k2(x,x')  is  the  second  iterated  kernel. 
Lonseth  [_!]  considered  the  example 


$ (x)  + / k (x,y)  $ <y)dy  = x2 , 0 £ x £ 1 , 


k(x,y)  = x(l-y)  if  x £ y,  k(x,y)  = y(l-x)  if  x £ y 


(5.10) 


for  which 


f (x)  = x2,  A = 1,  M - it"2,  L = 0,  5 = T'-- 1 -2-  , n = -~- 

x +1  x +1 


Kf(x)  = , K2f(x)  = -X~3|p-X  , ||  k(x,x')||x 


X * (1-x1) 


(5.11) 


The  crude  bounds  in  (5.1)  give 


0 < U(x')  - [x 


,2  _ 1]  < (15)_1/2  nx'  (1-x1)  , 0 < x-  <1 


(5.12) 


At  x'  = 1/2,  the  least  accurate  value,  this  shows  that 


0.21225  < $(1/2)  < 0.21819  . 


(5.13) 


The  bounds  (5.5),  with  $ proportional  to  f and  V zero,  give 


0.21704  < $(1/2)  < 0.21723  , 


(5.14) 


whereas  following  Lonseth's  approach  via  (5.4)  one  obtains  only 


0.21292  < $(1/2)  < 0.21740 


(5.15) 
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with  the  same  type  of  ♦ . Finally,  with  ♦ proportional  to  f and  ¥ proportional  to 
k,  the  more  sophisticated  bounds  (5.8)  yield  the  accurate  result 


0.217047  < $ ( 1/2)  <_  0.217048  . (5.16) 

The  integral  equation  (5.5)  can  be  solved  exactly  if  it  is  converted  into  a differential 
equation,  and  in  fact 

♦ (1/2)  = | sech(l/2)  - 2 = 0.2170472  . (5.17) 

6.  Bounds  on  Analytical  Approximations 

6.1.  The  Neumann  approximation 

Let  us  make  the  choice  of  trial  ve~tors 

*(x)  = 4>  (x)  , ¥(x,x’)  = V (x,x* ) , n,m  = 0,1,2...  (6.1) 

n m 

in  the  functionals  (4.1)— (4.3) , where 

tQ(x)  = 0 , 

♦ (x)  = {1  - AK  + xV  + ...  + (-AK)n  *}  f(x),  n > 1 , 
n — 

and 

¥0(x,x1)  = 0 , 

¥ (x,x* ) = {1  - AK  + A2K2  + ...  + (-Ak)1""1}  k(x,x' ) , m > 1 , 

m — 

= k(x.x')  - Ak  (x,x' ) + A2k  (x,x’>  — + (-A)m  1 k (x,x‘)  , 

2 j m 

the  functions  k_  ...  k being  the  iterated  kernels  (k,  = k)  . This  choice  leads  to  the 
2 m 1 

expressions 

J ( x ' ) = 0,  m + n - 0 , 

J(x')  = U - AK  + ...  + (-AK)m+n"1}  Kf  (x ' ) , m + n ^ 1 , 

S(x')  = (-AK)m+n  Kf(x')  (6.5) 

and 

C(x')  = | A | m+n  ||  k"  f ||  II  km+1<x,*'>llx  • <6-6> 

From  (4.8),  (4.9)  and  (2.17),  these  results  indicate  that  the  Neumann-type  approximate 
solution 

♦ <x')  - {1  - AK  + ...  + (-AK)“+n  + U-AK)m+n+1}  f(x’),  a < x'  < b , (6.7) 

is  in  error  by  not  more  them 


9- 


N 


m-tn+2 


(J  ♦ AM)  ( 1 - A 1 . ) 


l|Knf|!  I!  kml(*,x')||x 


The  relation  (2.6)  connecting  ij.(x)  and  i|-(x,y)  indicates  that  n = mil  ir.  a sensible- 
condition  to  impose,  but  the  maximum  error  in  (6.8)  holds  good  for  arbitrary  n and  m . 

6.2.  The  |N/N)  Tilde  approx i maul 

Consider  trial  vectors  of  the  form 

N-l 

*(x)  = l a Knf(x)  , (6.0) 

n=0 

N-l  N-l 

¥(x,x')  = J b Kn  k(x,x')  = T b k (x,x')  , (6.10) 

*•  n L‘  n n-fl 

n=0  n--0 

where  the  a and  h are  disposable  const  ants.  The  functional  J(x')  in  (4.1)  i r, 
n n 

itself  a bivariat ionol  approximation  to  X ^ { f ( x * ) - <(x')}  , and  so  it  is  not  un- 
reasonable to  choose  the  constants  a and  b to  make  J stationary  wits  respect  to 
variations  in  them.  Tf  this  is  done,  the  function  (f(x*)-  Xl(x')}  can  b identifi*  ! . 
the  (N/kJ  Pad'  approximant  to  4>(x')  in  powers  of  X,  constructed  directly  from  t:.- 
Neumann  Series.  Ch?sholm  [4]  has  suggested  this  sequence  of  apj  roxir.at  ions  to  * a:  nj. 
alternative  to  the  Fredholm  ones,  and  has  proved  a convergence  r«  suit  for  them.  The 

theory  of  the  present  paper  gives  bounds  on  them,  via  (4.8),  verifying  their  accur  :.cy  to 

, . 2M+1 

order  X 

If  we  use  the  notation 

f (x)  = Knf(x)  (6.11) 

n 


and  denote  by  the 

N 

array 

~\ 

(fi  4 

Af2) 

(f2  4 

•••  «N  + 

XfN+l» 

I 

(f2  4 

Af3) 

(f3  + XV 

(fN+l 

+ ' 

^ (6 

(fN  + 

XfN+l 

> <fN+l  4 

XW  ••• 

(f2N-l  4 Xf2N>  J 

then  the  functionals  J(x')»  S(x')  and  C(x')  can  be  v:r.i ‘.ten  as 


0 


J(x')  = - 


fNlX'> 


VX'> 


Vx*> 


(6.n) 


and 


C(X')  = ||  fiNf||  II  HNk(x,x')  ||  x , S(x')  = SiN  Kf(x')  , 


(6.14) 


where  the  operator  ft  is  specified  by 


X 


1 f1(x1) 

K f 2 (x * ) 


V*’1 


fNH(X,)  f2N(X'> 


v |an(x-)|  . 


(6. 1A) 


(cf.  Barnsley  (5)).  We  note  that  when  k(x,x')  is  a ('.eo-.-nerat.  hcrr.ol  of  order  the 

rows  of  tla  numcralor-deterni  n.u.t  in  (6.1I>)  an:  linearly  dependent,  and  ro  . ^ is  Uvn 
the  null-operator.  Acf  ordinal  the  [M/HI  bade  aj  . roni'nant  to  * is  exact  MJ . 

Bounds  on  other  types  of  Bade  appro/. ima.it  can  l,e  obtained  by  modifying  the-  pc.  rs 

of  K in  (6.9)  and  (6.10)  (cf.  6). 


6.3.  The  Fredholm  approxiiaat.io 


th 


One  choice  of  tri  al  vector''  whicti  leads  to  bound:,  on  ‘.lie  n Frr  al  o.m  ap;  : ■ ' - ' : * *'• 


*(>:)  - 

, ,n-l 

E i?.E,  + ...  + X E , 
_(' 1_  P-1 

f (x)  , 

n = 1,  2,  . . . ) 

d +ld  d . . . +)"d 

L 0 1 n _l 

J 

(6.16) 


1’(x)  “ 0 , 


where 


K - 

d 1 

- KE  _ , 

in  « l 

1U 

m 

n- 1 

m 

, . , . 1 

rod 

= l 

t d (-1)  , 

m = 1 

m 

r-1 

r m-r 

(h.  1 /) 
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and  t is  the  rth  trace  of  the  kernel  k(x,x')  . The  operators  E are  polynomials 
r m 

in  K#  and  KE  has  as  kernel  the  classical  Fredholm  minor;  the  numbers  d are  the 

IQ  IQ 

classical  Fredholm  determinants  • 

With  the  choice  (6.16),  f(x')  - Xj(x')  is  actually  the  nth  Fredholm  approxi- 

mation to  ♦(x')«  and  we  obtain  the  expressions 

J(x-)  = K4>(x’)  , 

C(x')  = ^ ——  ||(d  - KE  ) f !|  ||  Mx.x1)!!  (6.18) 

|d0+xd1+...+xndn| 

and 

XnKE  f(x') 

S(x')  = — . (6.19) 

d +Xd  +...+X  d 
0 1 n 

Similar  results  can  be  obtained  for  the  modified  Fredholm  approximation  should  the 
first  trace  t^  be  infinite. 

7.  Nonsyinmetric  Kernels 

If  k(x,y)  is  not  symmetric,  then  the  operators  K and  A = I + XK  are  not 
self-adjoint.  However,  if  A is  still  bounded  below  away  from  zero  in  the  sense  that 

II A 4> II  >_  y||  ♦Hi  Y > 0,  for  all  f e 3C  , (7.1) 

the  bivariational  bounds 

JCf,*)  - Y-1  II A 4—f  II  II  A*Y-g||  <<g,$>  £J(t,t)  + y”1  ||  A *-f  ||  ||  A*¥-g||  (7.2) 

* * * 

can  be  derived,  A = I + XK  being  the  adjoint  of  A . When  A = A , these  bounds  are 

coarser  than  those  in  (3.13).  Whenever  (7.1)  holds,  it  is  possible  to  infer  pointwise 
bivariational  bounds  on  4>(x')  by  taking  g(x)  = k(x',x),  although  the  special  role 
played  by  equation  (2.7)  specifying  the  reciprocal  kernel  is  lost.  In  (7.2),  V(x,x') 

would  not  now  approximate  the  reciprocal  kernel  of  k(x,x'),  but  rather  that  of  k(x',x)  . 
Details  are  given  elsewhere  [9] . 
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